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Abstract 

A coupled-channel model with ttN, pN and J/^N channels is developed to predict the 7r + N — > 
J/^ + N cross sections. The J/^-N interaction is parameterized in a form related to what has been 
predicted by the effective field theory approach and Lattice QCD. The other interactions within 
the model are constrained by the decay width of J/^f — > p + ir and the total cross section data of 
irN reactions. The calculated meson-baryon amplitudes are then used to predict the cross sections 
of the J/ty production on the deuteron target by including the contributions from the impulse term 
and the one-loop calculations of the final NN and J/*$>N re-scattering effects. Predictions of the 
dependence of the cross sections of 7r~ +p — > J/^ + n, 7 + c? — > J + n+p, and ir + d — > J/^ +p+p 
on the J/ty-N potentials are presented for experimental determinations of the J/^f-N interaction. 
Within the vector meson dominance model, we have also applied the constructed coupled-channel 
model to predict the 7 + p — > J/^ +p cross sections near the J/^ production threshold. 

PACS numbers: 25.20.Lj, 24.85. +p 
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I. INTRODUCTION 



The J/\& meson, as a cc bound state, can only interact with the nucleon through the 
gluon-exchange mechanism. By using the effective field theory method arid Lattice 
QCD jH[, it was found that the J/^-N interaction is attractive. Strong attraction was 
also foundjfl within a Pomeron-quark coupling model of cc-nucleus interactions. Thus it 
is possible that J/\l/ and a system of nucleons can form nuclear bound states with hidden 
charm. If such exotic nuclear states indeed exist and can be detected, we can get useful 
information for understanding the role of gluons in nuclei. 

The strength of the J/^-N interaction, parameterized as a Yukawa form vj/^n,j/^n = 
—a^-^-, deduced from Refs.[l|-[6| are rather different from each other. With p = 0.60 GeV, 
the calculated J/^-N scattering lengths range from -0.05 fm from using a = 0.06ofRef.fl to 
-8.83 fm from a = 0.60 of Ref.[6|. To understand quantitatively the role of gluons in hadron 
interactions, it is important to investigate how the J/ty-N interaction can be extracted from 
experiments to test these theoretical models. This is the objective of this work. Clearly, 
this is also a necessary step toward searching for possible nuclear bound states with hidden 
charm, as investigated in Ref. 0]. 

We first consider the it + N — > J/^> + N reaction. At energies above the J/^> production 
threshold, this reaction is strongly influenced by many possible reaction channels, such as 
pN, ttttN, irir-irN, nnKA channels, because of the unitarity condition. To proceed, we cast 
this complex reaction problem into a manageable coupled-channel model with only ttN, pN, 
and J/^/N channels. The interaction vj/^nj/^n is chosen to be of a form related to the 
Yukawa form of Refs.[l|-[6|. With the J/^-ir-p coupling constant determined^ by the par- 
tial decay width of J/\I/ — > irp, we can calculate the irN,pN — > J/tyN transition potentials 
v p n,j/vn and Vknj/vn by using the one-7r and one-p exchange mechanisms, respectively, as 
illustrated in Fig HI The other interactions, which do not connect with the J/^/N channel 
directly, will be treated as phenomenological complex potentials with their parameters con- 
strained by reproducing the total cross sections of irN reactions and ttN — > pN reactions in 
the considered energy region. We also constrain these potentials by reproducing the total 
cross section data of 7p — > p°p which can be related qualitatively to p°p —> p°p within the 
vector meson dominance model. To facilitate the experimental determination of the J/ty-N 
interaction, we will present results showing the dependence of the calculated cross sections 
of 7T~ + p — > J/ty + n on the potential wj/*at,j/*tv. 

With the meson-baryon amplitudes generated from the constructed coupled-channel 
model described above and the Pomeron-exchange model of 7 + N — > J/ty + N we had de- 
veloped in Ref.0], we then examine the dependence of the 7 + c? — > J +p + n cross sections 
on vj/q,N t j/^N. We also follow the suggestion of Ref. jsj to investigate ir + + d — > J + p + p 
reaction. Instead of using the factorization approximation used in Ref. 0, we have performed 
complete one- loop calculations to account for both the on- and off-shell J/^/N final state 
interactions. In addition, we also include the effects due to the final NN interaction by 
using the scattering t-matrix generated from the Bonn potential 

If we use the vector meson dominance hypothesis to convert the incoming photon into 
J/ty and p, we can predict the 7 + iV — > J '/^> + N cross sections within the constructed 
coupled-channel model. We will see that our predictions also depend on the J/^f-N potential 
vj/^n,j/^n- We will also present these results which could be tested in the forthcoming 



experiment [29 



In Section II, we present our coupled-channel model for J/ty production on the nu- 
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FIG. 1: Reaction mechanisms of 7r + N ->■ + iV(left) and p + iV -> J/V> + iV (right). 



cleon.The formula for calculating the meson-baryon reaction amplitudes are given in sec- 
tion III. In section IV, we give the formula for calculating the cross sections of the 
it + N — > J/fy + N and 7/7? + d — >■ J/ty + N + N reactions. Our results are presented 
in section V. A summary is given in section VI. 



II. COUPLED-CHANNEL MODEL FOR J/* PRODUCTION 

At energies near the J/^> production threshold, the ttN reactions involve many meson- 
baryon channels such as ttttN, ttttitN, ttttKA etc. Within the formulation of Ref.(9|, the 
scattering T-matrix for such reactions are defined by the following coupled- channel equations 

T a ,p(E) = V a j + ^ v a,-yGy(E)T^ t fi , (1) 

7 

where a,/3,7 denote the considered channels, G a (E) is the meson-baryon propagator of 
channel a, and v a ^ are the interaction potentials. 

For investigating the J/^-N interaction, we will treat ttN and pN channels explicitly 
since their interactions with the J/^/N channel can be calculated, as will be explained later. 
We thus set a, ,5,7 = J/^N,irN, pN, X in Eq.([T]), where X denotes collectively all other 
channels. By using the standard projection operator techniquejsl, [Toj], we can cast Eq.flT]) 
into 

T id (E) = Vij(E) + Vi,k{E)G k (E)T k j(E) , (2) 

k 

where k = J/^N,irN, pN. The energy- dependent interactions in Eq.([2]) are 

V^E) = «y +Y,%xG x (E)[l + t x , x (E)G x (E)}v X j , (3) 
x 

where the scattering amplitude t XtX in the subspace of the channel X is defined by 

t xx (E) = v x , x + v x , x G x {E)t xx (E) . (4) 

We now turn to defining the interaction potentials Vij in Eq.([3]). As noticed in Refs.0, 0], 
the decay width Tj/q, tP7T of J/^fr — > irp is significant. From the value of Tj/^ iP7T listed by 
the Particle Data Group (PDG), we can determine^ the coupling constant gj/^ t7Tp of the 
following interaction Lagrangian 



9j/v,-kp 



m 



(5) 
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We find[3] gj/^,n P = 0.032. By using Lj/^ np and the well known (TO] Lagrangian L n NN and 
LpNN, we can calculate one-7r (one-p) exchange transition potential of the p + N J/ty + N 
(tt + N — y J/fy + N) process, as illustrated in FigJT] Explicitly, their matrix elements can be 
written (with the normalizations to be specified in section III and omitting spin and isospin 
indices) as 



< P.J/ipiPN f \Vj/^N,TrN\PTT, f>N z > 



(27T) 



m N 



m N 



p7T 



E N (p Nf ) V E N (p Ni ) y/2E*{p*) 
1 



xu 



PN f 



9pNN ) Cfiuapq^Pj/TpZ (PJ/v)-^ 



mi 



m- 



4m 7v 



(6) 



< Pj/ip,PN f \Vj/ipN,pN\Pp,PNi > 



(2vr) 3 y/2Ej,^ m 

9j/ippw fnNN 



m N 



m N 



E N (p Nf ) V E N (p Ni ) ^/2E p (p p ) 



1 



e^P tl p £ u (Pp)Pj /lp e P (P 



x 



m 



2 U PN f 



u 



PN, J 



(7) 



where q = p^ — p^ f , the coupling constants g P NN = 6.20, k p = 1.825, f n NN = \/47r x 0.80 
are taken from a dynamical model Uj of irN scattering. All external particles of the matrix 
element of the above equations are on their mass shell with their momenta defined as p a = 
{E a {p a ),p a ) and E a {p a ) = [m 2 a + p 2 } 1 ' 2 . 

We assume that the mesons and baryons in the channel X do not contain charmed quarks. 
The coupling between X and J/^/N can then be neglected and we can set fj/vj/7v,x — 

0. It 

follows that for the interactions involving J/tyN, the second term of Eq.([n]) vanishes and we 
have the following simplification, 



VnN,J/^N 

v pN 

Vj/q,N,J/^N 



v ttN,J/^N , 
v pN,J/^!N , 
Vj/<f>N,J/^N 



The matrix elements for v^w jj^n and v pN)J /^, N have been given in Eqs. 
Vj/<5>N,j/yN, we follow Refs.[2j-l5j and use the following Yukawa form 



v j/vN,j/VN(r) = -a- 



-IM>r 



(8) 
(9) 
(10) 

and 0. For 



(11) 



To be consistent with the relativistic expression Eqs.(|6]) and ((7|), we assume that Eq.fflTT) is 
the non-relativistic limit of a one-scalar meson exchange amplitude in field theory. We then 
obtain 



<Pj/Tl>,PN s \vj/ii!N,J/ipN\Pj/y,PN i > 



(27T) 



m N 



m N 



1 



2E J/9 (pj /ip ) V E n(Pn s ) V E N (p Ni ) ^2Ejy{pj/v 
V 



x [gu^ip'j^e^p'j^)} 



2 U 'PN f U PN i 1 



(12) 
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where Vb = —&cmmj/^, and q = p^ i — p^ f . 

Since ivtv.x, v p n,x and txx(E) can not be calculated theoretically, we will determine 
Vij(E) with i,j = irN,pN in Eq.Q phenomenologically. For simplicity, we assume that 
they all have the following local form in the non-relativistic limit 

\ ;.,(/•) r'rH-:)J)Jr) ; i,j = nN,pN, (13) 

where fi,j{r) is unitless such as fij(r) = 1/(1 + e^ r ~ ro ^ t ) or e~ r2 l h2 . In the strong absorption 
(diffractive) modeljioj], the form of fi,i(r) for z = 1,2 is similar to the nucleon density. We 
will specify fij{r) later. 
We next define 

F M (g, E) = v°j(E) J e-^%(r)df . (14) 

Including appropriate covariant spin factors, the matrix elements corresponding to the form 
of Eq. (fT3~j) are taken to be 

< P , ^PNf\V 7 rN^N{E)\p^,p Ni > = J^y^pnj E nNtnN (q,E) up N . , (15) 

< i? P ,PN f \V pN , pN (E)\p p ,p Ni > = —--e v p {p' p )u ?Nf E pN , pN (q,E) u fN .s p>u (p p ) , (16) 

<^p,^iV / |V r pJV,7riv(^)|P 7 r, PNi > = / Q N 3 Pn £ »{Pp) E pN)7TN (q, E) up N [—J B ^Up N (17) 

[jl,7T J Win 777- 

We will determine the parameters v ? ■ and fij(f) phenomenologically by fitting the total cross 
sections of nN and nN — > pN, and the 7iV — > pN cross sections which can be calculated 
from the predicted pN — > pN amplitudes using the vector meson dominance hypothesis. 



III. CALCULATIONS OF REACTION AMPLITUDES 

In this work we follow the formulation of Ref. j9[ within which the scattering T-matrix is 
related to the S'-matrix by 

S U (E) = 6 U - (2n)t5(E f - Ei)6(P, - ?i)T u {E) , (18) 

where E a and P a are the total energy and momentum of the state a. The normalizations 
for the plane wave state \k > and bound state |$ a > are defined by 

< k\k' > = 8(k-k'), (19) 
< $ |$ 6 > = 5 ab . (20) 



A. Meson-Baryon reaction amplitudes 

In our calculations, we need the meson-baryon (MB) matrix elements in a fast moving 
frame. Following the instant form of relativistic quantum mechanics 12| , we write for the 
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two-particle M(pm) + B(p B ) M'(p M ) + B'{p B ) transition 



< PM m 'j M ^ m i M ',PB m 'j B ^ m r B \TM'B',MB(E)\p M m ™j s «V s > 



= J(P,k;p M ,PB)J(P ,k ;p M ,p B ) 
x <m' jM ,m' iM ;m' jB ,m' T Jt M 'B'MB(k\k]W)\m jAn m iM ;m jB1 m TB >, (21) 

where [rrij M , rrii M ] and [m Jfl , m TB ] are the z-components of the spin-isospin quantum numbers 
for M and B, respectively. The total momenta P, P ' and the energy W in the center of 
mass system are 

— * 

P = Pm + Pb , 
P ' = Pm + Pb , 

W = [E 2 -P 2 f 2 . (22) 

The relative momentum k and Jacobian J(P, k;pM,PB) i n Eq. ([2I|) are defined by the Lorentz 
boost transformation 

M M + Hq 

Tfpk-jt jM-1 |i/2 _ r E M (k)E B (k)H 1/2 

<J{^,k,Pm,Pb ) — \^f=, =r-r — F 7^ n F 7^~\JT l > l J4 J 



#o = E M (p M ) + E B (p B ) , (25) 



with 



M = \JH 2 -P* 

= E M {k) + E B (k) .) (26) 

In Eq. fl2ip . tM'S',MB(^) is the scattering operator in the center of mass frame. In the 
partial-wave representation, we can write 

t M 'B',MB(k', k; W) = Y \yi%?£w$ ') > tJ Js>,M'B',LS,MB( k 'i K W ) < ylfMB^iH 

JMj,TM T L'S',LS 

(27) 

where the angle-spin-isospin vector is defined by 

\yif, J M T B MT Ck)>= Y Y \ m jMi m iM'i m jB,m TB > 

3B m iM m,T B rni M 

[ ^ < J Mj\L S m L m s >< S m s \j M j B m hn m jB > 



3B 3M T B l M 

xl 

m L ,m s 

x < T M T \i M r B m iM m TB > Y L mL (k)] . 

(28) 



Here [jm, %m] and \j B , t b ] are the spin-isospin quantum numbers for M and B, respectively; 
< jm\jij2mim 2 > is the Clebsch-Gordon coefficient and Yt, mL {h) the spherical harmonic 
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function, as defined in Ref.[13(|. With the same partial- wave expansion Eq. (l28|) for the 
meson-baryon potential Vij(E), the Eqs.Q and (}2~Tj) lead to the following coupled- channel 
equation 

^l/S',m'b',ls,mb(^' \k.W) = V L ? s , M , B , LSMB {k',k.E) 

+ / k " 2dk " V L' T S',M'B',L"s",M"B"( k '^ k "- E ) G M»B" S " ,M" B" ,LS,M B ( k ^ kW ) 

L"S",M"B" J 

(29) 

where MB, M'B', hi" B" = ttN, pN, J/VN, and the propagator is 

1 



G MB (k,W) 



W - E M (k) - E B {k) + ze' 



We use the procedures developed in Ref. jl] to calculate the potential matrix elements 
^US' M'B 1 l"s" M"B"( k k" -E) from the matrix elements Eqs.fl6D-fTTTl). 

B. Deuteron wavefunction and NN amplitude 

For the calculations on the deuteron target, we need the deuteron bound state \^?p dt Mj > 



moving with a high momentum p^. Following Ref. [12|, it is defined by 



< pi, m Sl m T1 ;p 2 , m S2 m T2 \^^ Md > 

= 8{p d - p x -p 2 ) < p 1 ,m sl m T1 -,p 2 ,rn S2 m T2 \$p d7 M d >, (30) 

where 

< Pi,m Sl m Tl ;p 2 ,m S2 m T2 \^f di M d >= J(pd, «;Pi,P*2) < m Sl m Tl -,m S2 m T2 \x JdMd ' TdMTd (^) >(31) 
with 

\x mMT <*) >= E l^ dMTd (-) > • (32) 

Z=0,2 

Here = 1 is the deuteron spin and ui(k) is the usual deuteron radial wave function in its 
rest frame. The Jacobian J(p<i, n;pi,p 2 ) and the relative momentum k can be calculated by 
using the same Eqs.f l23|) -f l26|) with the replacement k — > k, pu — > Pi, Pb — > Pi, P ~^ Pd, and 
MB -»■ AT AT. 

We also need the NN amplitude 

< p^m' jv m' Tl )p 2 m' j2 ,m T2 \T NN , NN (E)\pxm h ^ > 
= J(P, k;p h p 2 )J(P' ,k';p[,p 2 ) 

X < m' jv m'^m' j2 ,m'^\t NN)NN {k ,k;W)\m jl ,m Tl ;m j2 ,m T2 > . (33) 



The above matrix element can be calculated from the same equations Eqs. (l2Tj) - (l28]) with 
the replacement of MB — > NN. We generated the NN partial-wave matrix elements 
^i!s'NN,LSNN(k' , k; W) from the Bonn potential [l^]. 
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IV. CALCULATIONS OF CROSS SECTIONS 



A. Cross sections of Meson- bar yon reactions 

We need to evaluate the cross sections for the reactions involving three meson-baryon 
(MB) channels with MB = nN,pN, J/^N. With the definitions Eqs.(fTg])-(12Tl). the differ- 
ential cross sections in the center of mass (CM) for the M(k) + B(—k) — > M'(k ) + B'(-k ) 
reaction can be written as 

da (An) 2 PM'B'(k')pMB(k) 



— V T \ 

dSl k 2 (2j M + l)(2j B + l) ^ 1 

I < m 'jM^ m iM^ m 'iB^ m 'r B \ t M'B',MB(k\k]W)\^^ > | 2 ] 



(34) 



where the matrix element of t j.f B' ,MB(k ; W) can be calculated from Eqs.( )27|) -( l28|) and the 
solution of Eq.( l29|) . The density of state is defined by 

kE M (k)E B (k) 

Pmb(k) = vr — . (35) 

The total cross sections a^ l N7rN of elastic nN — > nN and ct^n^n of nN — > pN can be 
calculated from using Eq.( )34|) . Eq.( )27|) and the solution of the coupled-channel equation 
Eq.®. 



B. Total cross sections of ttN reactions 

The total cross sections cr^ can be obtained from the irN elastic scattering amplitude 
using the optical theorem. With our definitions Eqs . (1 1 8 [) - f l2Tj) . we have 

a t ^ = -Im[f 7rN , 7rN (6 = 0)} (36) 

where the spin-averaged irN elastic scattering amplitude is 

Un,-kn(0) = < Q i m in'i m jNi m T N \UN,TrN(k\k;E)\0,m iw ;m jN ,m TN >] 

m m 

(37) 

and k = \k\ = \k'\ is the on-shell momentum, and cos9 = k • k. 

C. Cross sections of photo-production of J/^ and p 

We now note that with the vector meson dominance hypothesis, we can predict the cross 
section of 7 + p — > J/ty + p within the constructed coupled-channel model. This is done by 
writing the amplitude for 7(g) + p(— q) — > J/^(k ) + p(—k ) as 



< ^J/* ; m j Jv I^J/*P,7p(^''^ H/ )l A 7' m iiV > 
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= V < m'j /<s „m' jN \tj /9pt j /q , p (k' ,k; W)\m J/9 ,m jN > 5 m// ^,A 7 

+ < m 'j/y, m 'jjtj/*p,p p(k' ,k]W)\m p o,m jN > ydm^ , (38) 

where fj/q, = 11.2 and f p = 5.33. 

Obviously we can get the 7(g) + p{— q) — >■ p°(k') + p(—k') amplitude from Eq. (l3"8"|) by 
interchanging J/ty and p°. This will allow us to calculate the total cross section cr lp , p o p of 
IP ~~ >* P°P- Clearly, within this model based on the vector meson dominance hypothesis, 
&yp,p P is closely related to the total cross section of p°p — > p°p which can not be obtained 
experimentally, but can be essential in determining the parameters associated with our 
phenomenological potential V p n, p n(E). 




X 

N J/LfJ 



(c) J/ip N Re-scattering 

FIG. 2: Graphical representation of J/ip photo-production on the deuteron. Graphs(a), (b), and 
(c) correspond, respectively, to the impulse contribution, Eq. (fl3l) : NN re-scattering, Eq. (|50l) : and 
J/ipN re-scattering, Eq. (Tl9|) . 



D. Cross sections for 7 + ri — > J/^f + n + p 



In FigfJl we illustrate the mechanisms included in our calculations of the cross sections 
of J/^ production on the deuteron target. Since we are mainly interested in the J/^-N 
interaction (Figj2j(c)), we will examine how the predicted cross sections depend on the J/^fr- 
N relative momentum in the final J/^f + n + p state. This can be done most effectively by 
considering the following differential cross section in the CM frame of the 7(g) + d(—q) — > 
J/^{k) + n(p n ) + p(pp) reaction 

d 2 a A E^q)E d {q) f _^ E n (p n )Ej^(k)E p (p p ) M \p p \ 

dn p d\Kj/^\ 1 n ' \q\W J ^ ^E^-Kj/^Ej/^Kj/^) W 



2(2 J d + _ 

A 7 , m d mj/^,m n ,m p 
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I < kmjw, v P m p , Vn m n \T Imp + f J ^ N + f NN \q A 7 , $_ g - m<i > | 2 , (39) 

where q is the photon momentum chosen to be in the quantization ^-direction, $p d , md is 
the deuteron with momentum pd = —q, p n , p p and k are the momenta of the final proton, 
neutron and J/ip, respectively. Kj/^ = (\k,j/^\,Q Kj/4: ) is the momentum of J/W in the CM 

of the J/m-N subsystem. The amplitudes f Im P, f NN , and f J ^ N are calculated from the 
impulse ((a)), NN re-scattering ((b)), and J/^-N re-scattering ((c)) mechanisms shown in 

Figia 

In Eq.( l39|) . W is the invariant mass of 7c? and M = Ej/^,(Kj/^)+E n (Kj/^) is the invariant 
mass of the J/ip-N system. The magnitude of the outgoing proton momentum p p can be 
calculated fromM and W by 

\(W 2 -m 2 - Ml) 2 -AmlMl] 1 ' 2 . (40) 



2W 



The direction of p p is specified by fl p with respect to the incident photon momentum q. The 
other variables k for the outgoing J/\l/ and p n for the outgoing neutron in Eq. (l3"9"|) can then 
be calculated from k,j/^ and follows: 

~ k - *>» + k (sot - ■ (41) 

ft = - Xj/f + A (-^^y - • (42) 

With Eqs. (|40|) - (l42|) . all kinematic variables for calculating the integrand of Eq. (l39|) are com- 
pletely fixed for a given Q p and kj/^. Our task is to explore at what Q p , the calculated 
differential cross section ,„ d lz — , is most sensitive to the amplitude T J /^ N for the J/ty-N 

re-scattering mechanism (c) of FigJ21 

In the following three subsections, we give formula for evaluating the matrix elements 
Q £ rjnimp^ rpNN ^ an( j jij/vn _ ^y e evaluated the Clebsch-Gordon coefficients associated 
with the isospin quantum numbers in Eq. (l2Tj) and thus all isospin indices are suppressed 
and the amplitudes are on specific charged states specified explicitly as n for neutron and p 
for proton etc. 

1. Impulse Amplitude 

The impulse amplitude of Eq.(l39l(Fig]27a)) can be straightforwardly written as 

< k mj/^, p p m p , p n m n \f Imp \q A 7 , > 
= ^2 [ < ^ m J/1» Pp m p\ T J/i>P,~ip\q \i -Q-Pn m a >< -q-Pn m a , p n m n \$-$ mi > 

+ < k mj/Tp, p n m n \Tj/^n, in \q A 7 , -q-pp m a >< p p m p , -q-p p m a |$_ fmd > 

J (43) 

where < p p ;m p , p n ;m n \^ Pd -m d > has been defined by Eq. fl3Tj) (omitting isospin indices). 
The 7 + Ni — > J/tp + Nf amplitudes in the above expression are taken from our previous 
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work [3] 

< k mj/v, pN f rn Nf \Tj/y N}jN \q A 7 , p N . m Ni > 



m N 



m N 



^ yj2E j MPN f ) V Mm) ^JM 

x [um Nf (Piv>)e*(fc, m J/q ,)(M^(k,p Nf , q,p Ni ) + M^(k,p Nv q,p Ni )\e v {q, \)u mN% (p N .)) , 

(44) 

where .M^f and .M^ are the contributions from the Pomeron exchange and ir exchange 
mechanisms, respectively, as shown in the Figj3] 



Y vwv^ 

N 



2 S 

o C 

o c 

^4 



— J/y Y ww^-^-t->-J/|(/ 

Y 7T 



(a) 



(b) 



FIG. 3: Reaction mechanisms of 7 + A — > J/ip + A: (a) Pomeron-exchange, (b) pion-exchange. 



For the Pomeron exchange, the amplitude can be written as: 

M^ u (k,p Nf ,q,p Ni ) = Gp(s,t)7^r(t,g), 

with 

a P (t)-l 



[j-Y' 'expj— [a P (*)-l] 



(45) 



(46) 



T^(t,q) = il2y/^c 



2/ig 



^i(*)W^-M,(47) 



where t = (p^ — PN f ) 2 , s = (q — PNi) 2 , (*p(t) = «o + oi'pt with a = 1.25 and a' P = 1/sq = 
0.25GeV-\ $ = 1.1 GeV, a cm = e 2 /4vr = l/v^37, & = 0.84 GeV" 1 , / J/v , = 11.2, and 
Fx(t) = (AM 2 N - 2M)/(AM 2 N - t){\ - t/0.71{GeV 2 )) 2 . For the vr exchange, the amplitude 
Mt v is: 



k j U vn \ e 9J/<f,p K a fwNN ( A 
M^(k,p Nf ,q,p Ni ) = x 



2 \ 4 



(48) 

where the / p = 5.33, gj/v^ = 0.032, and A = 2000 MeV. 

All of the parameters specified above for evaluating Eqs. fj4*6|) - fj4*8|) were determined in 
Ref . 0] by fitting the total cross section data of 7 + p — > J/tp + p up to the invariant mass 
W = 300 GeV, as shown in FigH 
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FIG. 4: The total cross section of 7 + p — > J/tp + p. W is the invariant mass of the 7p system. 
The red solid curves are calculated from the model of Ref . 0] • 



2. J/ipN re-scattering amplitude 

The amplitude f J ^ N in Eq.([39D ( FigJ^c)) is : 

< kmj/^, p p m p , p n m n \f J/i)N \q A 7 , -q m d > 

E E 



3/„* 



d 6 k 



m a , m b =-1/2,1/2 m*^= 



-1,0,1 



X 



TV - E N (p p ) - Ej/^k*) - E N (-p p - k*) + te 

x < k* m} /v ,, p p m p \Tj/^ lp \q A 7 , k* +p p -qm a > 

x < kmj/^, p n m n \Tj/^j/^ n (Wj/^ n )\k* nCj/y, k + p n - k* m h > 

x <k*+p p -qm a , k + p n -k* m 6 |$_ 9 - Ad >] + [n o p] } (49) 



where Wj/y n = Ej/^(pj/^) + E N (p n ), < p m, p , ra'l^;^ > has been defined by Eq.flM]), 
and < kmj/Tp, p n Tn n \Tj/^ ntJ / 9n (Wj/ 9n )\k* m}^, p - k* m b > can be calculated using 
Eqs. (l2T]) - (l29|) . The first term in the bracket denotes the 7 + p — > J/tjj + p and J /ip + n — > 
J/ip + n, while the second term denotes 7 + n — > J/ijj + n and J/-0 + p — > J /ip + p. 



3. NN re-scattering amplitude 

The amplitude of Eq. (l39|) (FigJI^b)) can be written as: 

< kmj/^, p p m p , p n m n \f NN \q A 7 , -5 m d > 

rf 3 p* 

TV - Ej/^{k) - E N {p*) - E N {-k -p*)+ie 
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x ^2 < k mj/^i P* rn*\Tj/y N)lN \q A 7 , p m > 

m*, m, m' =-1/2, 1/2 

x[< p p m p , p n m n \T n p jn p(E np )\p* m*, p m! >< p m, p m'|$_g- md >] , 



(50) 



where p = p p + p n -p* , P = k+p* - q, E np = E N (p p ) + E N (p n ), < p;m, p '; m'|<% d;md > 
has been defined by Eq. fl3~TT) . and < p p m p , p n m n \T np!np (E n p)\p* m*, p m! > can be 
calculated by usin g E q.( l33|) and the NN — > NN partial wave amplitudes generated from 
the Bonn potential 3] • 



E. Cross section of tt + + d — > J + p + p 



By replacing Tj/^ Nr/N by Tj^ Ny7TN and changing notations appropriately, the formula 
presented in the previous subsection IV. B can be used to calculate the differential cross 
section -^-^ — , of 7r+ + d ->■ J/^ + p + p. 



V. RESULTS 

Our first task is to determine the parameters of the potentials of the coupled-channel 
model presented in section II. The parameters for the coupling potential V^j/^tv with 
i = irN, pN, J/^N have been specified there. We determine the parameters of Vij = 
Vij(E)fij(r) with i,j = TrN,pN, as already given in Eq. fTT3|) . by fitting the total cross 
section data of the nN reaction, a^ NlvN of nN elastic scattering, a n N, P N of rcN — > pN, 
and o" 7 p >p 0p of 7p — > p°p. We achieve this by choosing 

^ = 1 + e (r-c)A ( 51 ) 

where c = 0.8 fm, t = 0.4 fm for all i, j = ttN, pN. The energy dependence of the total cross 
sections in the energy region near the J/\l/ production threshold (4 GeV < W < 6 GeV) 
can be fitted qualitatively by setting 

vUE) = (-z)[A + B(E-E,f] (52) 

Our fits to o l °tp are shown in Fig]5j The small isospin dependence in the data is neglected 
in our fits. In Figj6j we see that the fits to the total cross section data of 7r ± + p — > p ± + p 
(left) and 7 + p — > p° + p (right) are also reasonable. The resulting parameters A, B, and 
Eq are listed in Table [D 

Note that the forms Eqs.( l5ip -f l5"2"|) are purely phenomenological, and our fits in FigsJMni 
are only qualitatively. But they are sufficient for estimating the order of magnitudes of 
the J/ty production cross sections. For a more quantitative calculation, we clearly need to 
improve this phenomenological aspect of our model. 

With the fits shown in FigsJMSl all of the parameters for our calculations are completely 
fixed. In the next few subsections, we present our predictions for future experimental deter- 
minations of the J/ty-N interaction. 
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TABLE I: The parameters for v®j(E) of Eq. (|52p . E is the total energy in the enter of mass system. 

~~ i j A (GeV) B (GeV- 1 ) E (GeV) 
ttN ttN 0.43 1.18 4.0 

ttN pN 0.04 0.09 4.0 

pN pN 0.45 0.45 4.0 
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FIG. 5: The fits to the data of the total cross sections a to + of ir^p reactions, and the total elastic 



ir^p — > -K^p cross section a el + + . The data are from PDG 241 



n 



The simplest experiment to determine the J/^-N interaction within our coupled-channel 
model is to measure the cross sections of ir~ + p — > J /\& + n reaction. We first observe 
that the coupled- channel effects can change drastically the magnitudes and the energy- 
dependence of the predicted cross sections. This is illustrated in the left side of FigJTJ In 
the absence of coupled- channel effects, the results from the Born approximation (setting 



■ J/^N,ttN 



J/^N,ttN, 



are the dotted blue curve. When the interactions associated with 



the J/^/N channel, V^j/^ with i = irN, pN, J/tyN, are included in solving the coupled- 
channel equation, we obtain the dashed red curve which is suppressed at high W in contrast 
to the raising behavior of the dotted blue curve from the Born approximation. When the 
complex potentials Vij with i,j = 7iN,pN are also included in our full calculations, we 
obtain the solid black curve. We see that these complex potentials can reduce drastically 
the magnitudes of the cross sections. This is due to the fact that most of the incident pions 
are absorbed before the J/^> production takes place. Within our model, this absorption 
effect is due to the very large imaginary part of V^n^n- 

The coupled- channel effects on J/ty + p — > J/ty + p are shown in the right side of FigJTl 
We note that the J/^-N interaction, as defined in Eq. fllip . is real and therefore does not 
have strong absorption effects on J/ty-N scattering. By comparing the differences between 
the blue dotted curve of the Born approximation (Tj/^nj/^n = Vj/^n,j/^n) and the red 
dashed curve, we see that the coupled- channel effects due to V^j/-^ with i = nN,pN, J/^/N 
is to increase the cross section at all W. The coupled- channel effects due to Vjj with 
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FIG. 6: The fits to the data of the total cross sections cr n ± p ^. p ± p of ir^p — > p^p (left) and a. 



for 7p — > p°p (right). The data of cr- K ± p ^ tp ± p is from 25], and the data of o" 7pH>p o p is from 261428], 



i,j = irN,pN are very small, as can be seen by comparing the red dashed curve and the 
black solid curve in the right side of FigfTl This is due to the fact that the J/ty-n-p coupling 
constant gj/^ pn = 0.032, calculated from the partial decay width of J/^/ — > irp, is very small 
in the calculation of Eqs.0H])-([Z]) for the one-meson-exchange matrix elements of f p v,j/* and 




FIG. 7: The cross sections of 7r~ +p -> J/^ + n (left) and J/^f + p ->• J/^ +p (right). Black 
Solid: full calculation, red dashed: Coupled-channel, but setting V it j/^ = for i,j = 7rN,pN, blue 
dotted : Born calculation Tj = Vij. W is the total energy in the center of mass frame. 

We next examine the dependence of the ir~+p — > J +n cross sections on the strength a 
of the J/^-N potential vj/^n,j/^n- The results from using a = 0.2, 0.09, 0.06 are compared 
in the left side of FigJB] The predicted total cross sections at peak positions are about 1.5 
nb. Such small cross sections are due to the strong absorption of the incident pion and that 
the irN J/^N transition potential V n N,j/^N is weak, as discussed above. The differences 
between three results are significant only in the region near W ~ 4.25 GeV. On the other 
hand, the predicted J/^/N — > J/^/N cross sections are more sensitive to a, as seen in the 
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right side of FigJHJ This suggests that J/ty-N interaction can be more easily determined in 
the reactions of J/^f production on the deuteron target in the kinematic region where the 
cross sections are sensitive to the J/^N — > J/^/N re-scattering mechanism (Figj2(c)). This 
is what we will examine in the next two subsections. 




(a) W(GeV) (b) W(GeV) 



FIG. 8: The total cross sections of tt + p — > J/^f + n (left) and J ftp -\-p — > J/^> +p (right). W is 
the total energy in the center of mass frame. The black solid, blue dotted, and red dashed curves 
are calculated by using the J/V-N potential Eq.jTT]) with y, = 0.6 GeV and a = 0.20,0.09,0.06, 
respectively. 



B. 7 + d ->• J/* + n + p 

We only consider the energy region close to the J/ty production threshold. The results 
presented below and in the next subsection are calculated at the energy 100 MeV above the 
J/\& production. 

As discussed in section IV.D, our main task is to use Eq. (l39|) to examine the dependence 
of the differential cross section dn jjfj~^ °f this process on the J/ty-N potential Vj/vpjv.JW- 
Obviously, we need to identify the region of the outgoing proton angle 9 P with respect to 
the incident photon where the J/^-N re-scattering mechanism ((c) of FigfJ]) dominates. 
We find that this is in the region close to 9 P = 0, as can be seen in the left side of FigJHJ 
We see that in the low J/\& momentum, kj/<^ < about 200 MeV, the contribution from the 
J /ty re-scattering term (blue dotted curve) is much larger than that from the other two 
mechanisms. On the other hand, the impulse term dominates at 8 P = 180° as seen in the 
right side of Figj9j Accordingly, three J/^-N potentials with the strengths a = 0.2, 0.09, 
and 0.06 obtained by the effective field theory approach and LQCD can be easily tested at 
9 P = 0, but less possible at 9 P = 180°. This is shown in FigJTUJ 

The results presented above clearly indicate that it is highly desirable to have data in 
the region close to 9 p = in future experiments in order to be able to determine the J / ty-N 
interaction. 
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FIG. 9: The differential cross section da/ (dflpdnj/^,) of j+d — > J /^ +p+n vs the momentum Kj/q, 
of J/vl/ in the center of mass frame of the J/^+n sub-system. The results are from coupled-channel 
calculations using the J/^-N potential, Eq. (|lip . with p = 0.6 GeV and a = 0.20. 9 P in (a) and (b) 
are the angle between the incoming photon and the outgoing proton. The red dashed, blue dotted 
and pink dashed-dotted lines are the contributions from the amplitudes of the Impulse term, J/^N 
re-scattering, and NN re-scattering, respectively. The black solid lines are the coherent sum of 
these three amplitudes. 




FIG. 10: The differential cross section da /{drtpdnj/^,) of 7 + d — > J/ty + p + n vs the momentum 
of J/^f in the center of mass frame of the J/^f + n sub-system. 6 P is the angle between the 
incoming photon and the outgoing proton. The black solid, blue dotted, and red dashed-dotted 
lines are from the coupled-channel calculations using the J/^fN — > J/^N potential, Eq. (|11|) . with 
H = 0.6 GeV and a = 0.20, 0.09, 0.06, respectively. 



C. 7T+ + d ->• J /iff +p + p 



For the calculated dn | of the n + d — > J/ty + p + p reaction, we find that the impulse 

term dominates at all angles. Furthermore, the contributions from the J /^>-N re-scattering 
(Fig E(c)) are weaker than those of the NN re-scattering (FigJ2(b)). Nevertheless, the data 
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from experiments on this reaction can be useful to determine the J/ty-N interaction. This 
is illustrated in FigJTTI where the results from a = 0.20, 0.09, and 0.06 for 8 P = and 180° 
are compared. 




k j/v (MeV) k j/v (MeV) 



FIG. 11: The differential cross section da/(d£l p dKj/q,) of n + + d — > J/^ +p + pvs the momentum 
Kj/q, of J/\I / in the center of mass frame of the J/^f + p sub-system. The black solid, blue dotted, 
and red dashed-dotted lines are from the coupled-channel calculations using the J/^-N potential, 
Eq. tjlip . with \i = 0.6 GeV and a = 0.20, 0.09, 0.06. 6 p is the angle between the incoming pion and 
the outgoing proton. 



D. 7p — > J/*S>p in the coupled-channel model 

With the vector meson dominance hypothesis, we can predict the cross section of 7 + p — > 
J/ty +p within the constructed coupled- channel model. This is done by using the amplitude 
given in Eq. (1381) . For three coupled-channel amplitudes calculated with a = 0.2, 0.09, and 
0.06 for Vj/^at j/^at, our results are compared with that of the Pomeron-exchange model 
(green dot-dashed curve) in the left of FigJT2l 

We note that the coupled- channel model results are expected to be valid only at energies 
near J/\I/ production threshold where the data are rather uncertain. On the other hand, the 
use of Pomeron-exchange model at very low energies is also questionable. Our results for 
a = 0.06 are close to the three data points near 4.5 GeV. At W > about 5 GeV, we clearly 
need to improve the model. The data from the forthcoming experiment (29| at JLab will be 
useful to clarify the situation. 

In the right side of FigJT2| we show that the coupled-channel effects can increase the 
cross section significantly. Our results here as well as that shown in FigJT] suggest that any 
attempt to determine the J/ty-N interaction must include the coupled-channel effects, as 
required by the unitarity condition. 

VI. SUMMARY 

We have developed a coupled- channel model with ttN, pN and J/^fN channels to pre- 
dict the J/ty production on the nucleon. The J/^-N interaction potential Vj/^j/^tv is 
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FIG. 12: The cross section of 7 + iV — > J/ty + N reaction. W is the total energy in the center 
of mass frame. In the left side, the black solid, blue dotted, red dashed lines are calculated by 
using couple channel model with \i = 0.6 GeV and a = 0.20, 0.09, 0.06 and the green dashed-dotted 
line is from Pomeron-exchange. In the right side, the red dashed and blue dashed-dotted lines are 
from the full calculation and the Born approximation calculation Tj » = Vij. The parameters of 



v j/^n,j/^n of the coupled-channel model are p, = 0.6 GeV and a = 0.06. The data is from |l5l-ll7l]. 



taken from the calculations using the effective field theory method and Lattice QCD. The 
J/^N — > 7iN,pN transition potentials Vj/^N,nN and Vj/^n, p n are calculated from one-p- 
exchange and one-7r-exchange mechanisms with the parameters determined by the decay 
width of J/'ty — > up and the previously determined uNN and pNN coupling constants. 
The interactions Vij with i, j = ttN, pN are treated as phenomenological complex potentials 
with their parameters constrained by the fits to the data of total cross sections of the ttN 
reactions, irN — > irN, irN — > pN, and 7£> — > p°p. 

The calculated meson-baryon amplitudes are then used to predict the cross sections of 
7 + d J/^ + N + N and n + d J/ty + N + N reactions. The calculations on the 
deuteron target involve the contributions from the impulse term and the final NN and 
J/^N scattering. We have identified the kinematic region where the J/^-N potentials can 
be distinguished in 7r~ + p — > J/^ + n, 7 + c? — > J/fy + n + p and 7r + + <i — >■ J/^> + p + p 
reactions. Predictions of the dependence of the cross sections of these reactions on the J/^>- 
N potentials are presented. Our results shown in FigsJHJ [TOj and [11] can be used to facilitate 
the experimental determinations of J/^-N interactions. 

Within the vector meson dominance model, we have also applied the constructed coupled- 
channel model to predict the 7 + p — > J/^ + p cross sections. Our results near the J/^f 
production threshold, as shown in Figfl2| are very different from what can be calculated from 
the conventional Pomeron-exchange model which is mainly constrained by the 7p —> J/^+p 
cross section at high energies. It will be interesting to distinguish these two models in the 
forthcoming experiment at Jefferson Laboratory [29]. 
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